NASA  TR  R-8 


"  "s. 


/  “' "  .i 


TECHiNieAL  IMPORT 

V  :'.R-^  V^;  :'' 


Reproduced  From 
Best  Available  Copy 


By  FRED  W.  MAtTING,  DEljyN  ^  CHAPMAN/ JACk^  R 

/'  J  ’  '  ' vilia^ANDBiw.  d^THOMlAS'  ■,  ',  ■  - 


'  >/^  - 


V'  ,  (- 


.  .IVV 


yXQ,.^  i  C 


E«STl^lBirnOHSTAtE^^A 

A^prd^^dfoVpfiijild Release  o 
^  Distribljtton  ynllmited  i  >x  ^ 


196K,  < 


j-  ■'■  /-/- 


20011)05  069 


For  jmleliy  the^  S^erlntendeni  of  ItaeiUneAt^  U.S.  ^^enunent  Print^  Office,  Washloston  2S,  p.C.  Teorly  ^t^riplioii,  $15;  forei^,  |19 
<  /^  single copy  priw  vidoa iicww^iig  v  Prl^  ^ ^  ^\ 


TECHNICAL  REPORT  R-82 


TURBULENT  SKIN  FRICTION  AT  HIGH  MACH  NUMBERS  AND 
REYNOLDS  NUMBERS  IN  AIR  AND  HELIUM 

By  FRED  W.  MATTING,  DEAN  R.  CHAPMAN,  JACK  R.  NYHOLM,  and  ANDREW  G.  THOMAS 

Ames  Research  Center 
Moffett  Field,  Calif. 


TECHNICAL  REPORT  R-82 


TURBULENT  SKIN  FRICTION  AT  HIGH  MACH  NUMBERS  AND  REYNOLDS 

NUMBERS  IN  AIR  AND  HELIUM 

By  Fred  W.  Matting,  Dean  R.  Chapman,  Jack  R.  Nyholm,  and  Andrew  G.  Thomas 


SUMMARY 

Experimental  measurements  of  local  skin  friction 
in  the  turbulent  boundary  layer  have  been  made  over 
a  Mach  number  range  from  0.2  to  9.9  and  a  Reynolds 
number  variation  from  2X10^  to  100X10^.  Flows 
were  two-dimensional  over  a  smooth  flat  surface  with 
zero  pressure  gradient  and  with  adiabatic  conditions 
at  the  wall.  Air  was  used  as  the  working  fluid  up 
to  a  Mach  number  of  Jf. 2  while  helium  was  used  from 
a  Mach  number  of  4.2  to  a  Mach  number  of  9.9. 
Boundary-layer  profiles" at  selected  Reynolds  numbers 
were  taken  at  all  supersonic  Mach  numbers  at  which 
tests  were  run. 

An  equivalence  parameter  for  comparing  turbulent 
boundary  layers  in  different  fluids  has  been  used  and 
experimentally  verified.  By  means  of  this  parameter^ 
experimental  results  obtained  with  helium  have  been 
converted  to  equivalent  air  results. 

Several  methods  of  calculation  of  skin  friction  in 
the  compressible  turbulent  boundary  layer  have  been 
compared  with  experimental  results.  Of  the  methods, 
the  intermediate  enthalpy  (T')  method  appears  to 
give  answers  closest  to  the  experimental  results  over 
the  wide  ranges  of  Mach  and  Reynolds  numbers  used. 

INTRODUCTION 

The  history  of  research  in  the  turbulent  bound¬ 
ary  layer  can  probably  be  characterized  fairly  by 
the  statement  that  a  tremendous  quantity  of 
measurements  have  been  made,  but  there  is,  as 
yet,  no  complete  theoretical  formulation  of  the 
problem.  The  lack  of  a  satisfactory  theoretical 
treatment  has,  of  course,  created  the  need  for  the 
large  volume  of  measurements  that  have  been  and 
are  being  made. 

Formerly  the  principal  regime  of  interest  in 
turbulent  boundary  layer  research  was  in  low 


Mach  number  (i.e.,  incompressible)  flow.  The 
older  measurements  have  been  well  summarized 
by  Schoenherr  (ref.  1),  while  later  incompressible 
experiments  have  been  surveyed  by  Locke  (ref.  2) 
and  Hughes  (ref.  3).  The  most  recent  definitive 
measurements  in  low  Mach  number  turbulent 
boundary  layers  were  made  by  Smith  and  Walker 
(ref.  4)  who  employed  direct  force  measuring 
techniques.  It  can  be  stated  that  the  great  mass 
of  measurements  (particularly  skin  friction)  in 
low-speed  turbulent  boundary  layers  did  fill  a 
real  need  in  supplying  design  data  in  the  fields 
of  naval  architecture,  hydraulic  machinery,  and 
subsonic  aircraft. 

As  power  plants  capable  of  propelling  aircraft 
at  supersonic  speeds  became  available,  the  field 
of  interest  naturally  shifted  to  turbulent  boundary 
layers  in  supersonic  and  hypersonic  flows.  Com¬ 
pressibility  effects  in  the  high-speed  boundary 
layers  introduced  additional  variables  and  param¬ 
eters  to  the  problem,  thus  increasing  the  number 
of  boundary  conditions  for  which  measurements 
were  desired.  At  the  same  time  the  insertion  of 
new  variables  and  parameters  has  made  the  theo¬ 
retical  problem  even  more  formidable.  A  number 
of  measurements  of  compressible  turbulent  skin 
friction  made  in  the  last  decade  have  been  re¬ 
ported  in  references  5  to  20,  inclusive. 

Almost  all  of  the  above  referenced  measure¬ 
ments  were  taken  over  a  range  of  Reynolds  num¬ 
bers  and  Mach  numbers  which  is  small  compared 
to  that  encompassed  by  current  vehicles.  The 
present  program  of  research  was  initiated  at  the 
Ames  Research  Center  in  order  to  secure  data  for 
turbulent  skin  friction  at  high  Mach  and  high 
Reynolds  numbers  where  no  previous  data  have 
existed.  These  data  were  obtained  by  direct 
force  measurements  for  adiabatic  wall  conditions. 
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The  program  included  the  construction  of  a  special 


apparatus,  a  1-  by  10-inch  boundary-layer  chan¬ 
nel,  with  provisions  for  using  both  air  and  lielium 
as  test  media.  In  connection  with  using  helium. 

it  was 

necessary  to  determine  the  equivalence  of 

air  and  helium  in  the  turbulent  boundary  layer. 

This  is 

covered  in  some  detail  in  later  sections. 

SYMBOLS 

A 

area 

a 

speed  of  sound 

Of 

coefficient  of  local  skin  friction,  ,r  2 

Gn 

coefficient  of  local  skin  friction  for  incom- 

prcssible  flow 

Cp 

pressure  coefficient, 

Cp 

specific  heat  at  constant  pressure 

€p 

specific  heat  at  constant  volume 

h 

enthalpy  per  unit  mass 

k 

coefficient  of  thermal  conductivity 

L 

length  from  nozzle  tliroat  to  center  of 

skin-friction  element 

1 

characteristic  reference  length 

M 

Mach  number 

Ma 

equivalent  air  Mach  number  (see  eq,  (9)) 

V 

static  pressure 

Pr 

Prandtl  number,  ^ 

fuib 

c  c 

turbulent  Prandtl  number, 

Q 

convective  heat- transfer  rate  per  unit  area 

in  the  boundary  layer  in  the  y  direction 

Qturb 

turbulent  heat- transfer  rate  per  unit  area. 

in  the  y  direction 

dynamic  pressure,  ;  or  {\j2)pJJJ 

R 

gas  constant 

Rcl 

Ke^molds  number  based  on  L,  ^ 

Mm 

Rci 

Reynolds  number  based  on  1,  ^ 

Rex 

Reynolds  number  at  skin-friction  element 

measured  from  transition  point, 

U 00 Poo  (-^  ^tr) _ U 

Mco  Mm 

RCy) 

Reynolds  number  based  on  wall  properties. 

u^pj  /  u^p„x  A 

M«.  \  Mi^  / 

T 

static  temperature,  °R 

t 

static  temperature,  °F 

v,r 


Vr 


Ucx> 


O'lr 

V 

y* 

y 

A 

d 

€ 


Vr 

e 


X 

M 

V 

i 

r 

turl} 


oo 

t 


turh 

tr 

tv 

aw 

0 


velocity  components  in  the  x  and  y  direc¬ 
tions,  respectively 

friction  velocity, 

\  Ptr 

dimensionless  velocity,  — 

tij 


frec-strcam  velocity 

coordinate  parallel  to  free-stream  flow 
direction  and  Icngtli  from  transition 
point  to  center  of  skin-friction  element 
length  from  no^^zle  throat  to  transition 
point 

coordinate  perpendicular  to  the  wall 
dimensionless  y  coordinate, 


Q 

ratio  of  specific  heats,  —■ 

increment  of  a  quantity  (e.g.,  AM) 
boundary-layer  thickness 

eddy  viscosity, 

dimensionless  coordinate  perpendicular  to 
the  wall,  j 

.  .  Ta,r-T^ 

recovery  factor,  -j,  ^  j  ~ 


b  o  II  n  d  a  r y  -1  a  y  e  r  momentum 


eddy  difi'usivity, 
viscosity 


thickness, 


kinematic  viscosity,  ^ 

dimensionless  coordinate  parallel  to  frcc- 
stream  flow  direction,  j 


mass  density 

shear  in  the  boundary  layer 
turbuh'Tit  shear  in  the  boundary  layer 


SUBSCRIPTS 

free-stream  conditions  ai  element  location 
local  total  conditions,  for  example, 

turbulent 
transition  point 
wall  conditions 
adiabatic  wall  conditions 
gas  other  than  air 
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SUPERSCRIPTS 

'  reference  conditions  with  the  T'  method,  or 
fluctuating  quantities  in  turbulent  flow 
—  time  mean  quantities  in  turbulent  flow 
~  dimensionless  time  mean  quantities  normal¬ 
ized  with  respect  to  the  corresponding  free- 

stream  quantity,  for  example, 

CONSTANT-PRESSURE  TURBULENT  BOUNDARY 
LAYERS  IN  DIFFERENT  CASES 

ANALYSIS  OF  EQUIVALENCE 

The  purpose  of  the  analysis  which  follows  is  to 
develop  the  conditions  of  equivalence  for  turbu¬ 
lent  boundary  layers  in  different  gases  with  con¬ 
tinuum  flow.  To  do  this  the  equations  of  motion 
and  energy  will  be  expressed  in  dimensionless 
form  and  the  resulting  dimensionless  parameters 
and  boundary  conditions  studied. 

The  equations  of  the  constant  pressure,  steady 
state,  compressible  turbulent  boundary  layer  fur¬ 
nish  a  convenient  starting  point  when  put  in  the 
following  form  (ref.  21): 


Continuity: 

A(piz+7^)+l(pij+77)=o  (1) 

x-momentum : 

(ju+7^)  If +(p«+^)  1^=1 

(2) 

Energy : 

Equations  (1),  (2),  and  (3)  can  be  written  in 
dimensionless  form  if  all  quantities  are  normalized 
with  respect  to  free-stream  quantities  (e.g., 
p  — - pj pea }  ^ — u/U^co^  etc.). 

Continuity : 


x-momentum : 


/ —  ,  p'u'  \  ^  f 


where 


P»4  ^ 


p'  \  du 


~ _ ^ 


Energy: 

/ —  ,  p'u'  \  dh  ,  I  p'v'  \ 


brj 


Rei  dr)  L\ 


JI  \dh 


PvtuTb  PrJ  dr]_ 

These  dimensionless  equations  are  subject  to  the 
following  boundary  conditions  on  velocity  and 
enthalpy : 

o)=o  X({,  0)=!*, 


«>)  =  1  h{^,  co)=l 

Two  turbulent  boundary  layers  will  be  equiva¬ 
lent  (dynamically  similar)  when  the  dimensionless 
differential  equations  and  dimensionless  boundary 
conditions  for  the  two  boundary  layers  are  identi¬ 
cal.  Examination  of  equations  (4),  (5),  and  (6) 
shows  that  equivalence  can  be  expected  to  exist 
between  turbulent  boundary  layers  in  different 
gases  when  the  following  dimensionless  parameters 
and  variables  are  matched: 


Eei 

UJ/k 

hw 


Controllable  parameters,  required  for 
matching  certain  coefficients  of  the  differ¬ 
ential  equations  and  the  boundary  condi¬ 
tions  on  enthalpy. 


p(h)  State  variables,  required  for  matching  the 
Pr(fi)  thermodynamic  state  and  the  laminar 
JiCh)  J transport  properties. 


Turbulence  variables,  required  for 
matching  the  eddying  motion  rep¬ 
resented  by  the  fluctuation  terms 
in  the  differential  equations. 


P^  turb{.^jV^ 


As  listed  above,  the  parameters  and  variables  are 
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divided  into  three  groups:  the  controllable  param¬ 
eters  which  can  be  matched  by  control  of  test 
conditions;  the  state  variables  which  may  or  may 
not  be  matched,  depending  on  the  particular 
thermodynamic  and  transport  properties  of  a 
given  gas;  and  tlie  turbulence  variables,  the 
matching  of  which  must  be  determined  by  experi¬ 
ment.  Experimental  results  are  needed  since 
there  are  available  no  exact  analytical  expressions 
for  the  turbulence  variables  which  apply  through 
the  entire  boundary  layer.  (As  a  consequence 
there  are  no  analytical  solutions  of  the  system 
of  equations  (4),  (5),  and  (6).)  Although  the 
turbulence  terms  can  be  measured  from  point  to 
point  in  comparing  two  boundaiy  layers,  it  is 
probably  more  convenient  simply  to  measure  the 
over-all  effects  of  these  variables  on  more  readily 
measured  quantities  such  as  skin  friction  and 
boundary-layer  profiles;  the  latter  method  is 
used  herein. 

In  setting  up  equivalent  turbulent  boundary 
layers  in  two  different  gases,  some  of  the  para¬ 
meters  and  variables  can  be  easily  matched 
while  the  matching  of  otlicrs  presents  difficulties. 
The  controllable  parameters  can  be  exactly 
matched  between  two  turbulent  boundaiy  layers. 
The  three  controllable  parameters,  UJIh^, 

hy,,  can  be  thought  of  as  three  independent  dimen¬ 
sionless  energy  ratios  obtained  from  six  dimen¬ 
sional  boundary  values,  hrr,  1. 

The  Reynolds  number,  Bci,  represents  the  ratio 
of  inertial  to  viscous  energies  in  the  free  stream; 
UJIh„  is  double  the  ratio  of  kinetic  energy  to 

enthalpy  in  the  free  stream,  while  hy,  represents 
the  ratio  of  wall  enthalpy  to  free-stream  enthalpy. 
Even  without  recourse  to  the  differential  equa¬ 
tions,  it  is  clear  that  dynamic  similarity  between 
two  turbulent  boundary  layers  would  not  be 
possible  unless  these  energy  ratios  were  matched. 

The  matching  of  the  state  quantities,  p{h), 
PrQh),  p(h),  depends  on  the  gases  being  compared. 

a.  The  quantity,  p{h),  will  be  well  matched 
when  the  gases  being  compared  do  not  deviate 
greatly  from  perfect  gases.  For  thermally  perfect 
gases  at  constant  pressure : 


P(h) 


1 

dh 

_  1  J 

0  c\(h) 

m 

(7) 


If  the  gases  are  also  calorically  perfect: 


p(h)=^f 

h 


(8) 


b.  The  Prandtl  number  is  nearly  constant  and 
is  approximately  the  same  for  all  gases  in  the 

absence  of  imperfect  gas  effects;  hence  PrQi)  will 
generally  be  well  matched. 

c.  Tlie  dimensionless  viscosity,  p{h)^  is  a  param¬ 
eter  which,  in  general,  may  not  ])e  perfectly 
matched  between  two  difl'erent  gases.  However, 
tins  match  can  be  reasonably  well  approximated 
by  the  selection  of  gases  and  temperature  levels. 

There  is  no  effective  method  of  predetermining 
tliat  the  individual  turbulence  variables  can 
actually  be  matched  at  all  corresponding  points 
in  two  boundary  layers.  However,  it  is  known 
tliat  if  true  dynamic  similarity  exists,  the  dimen¬ 
sionless  differential  equations  (with  boundary 
conditions)  will  be  identical  for  two  boundary 
layers.  As  a  consequence,  the  coefficients  of 
local  skin  friction,  Cf=Ttpl{Vi)p^TU,  and  the 
dimensionless  velocity  profiles,  ft(|,  ??),  must  be 

the  same  for  the  two  boundary  layers.  ^The 

riglit-hand  side  of  cq.  (5)  may  be  written  as 

i  clclonnin- 

ing  that  the  C/s  and  the  dimensionless  velocity 
profiles  arc  identical,  it  can  be  concluded  that  a 
necessary  condition  for  dynamic  similarity  has 
been  fulfilled.  In  the  strict  sense  this  is  not 
sufficient  to  prove  dynamic  similarity,  but  one  is 
encouraged  to  believe  that  dynamic  similarity 
exists  in  this  case  for  all  practical  purposes. 

The  analysis  above  has  been  concerned  only 
with  turbulent  boundary  layers,  but  it  is  of 
interest  to  note  that  for  dynamic  similarity 
betweeii  laminar  boundary  layers  (at  zero  pres¬ 
sure  gradient)  in  different  gases,  the  quantises 
that  must  be  matclied  are  the  controllable  param¬ 
eters  and  tlie  state  variables  as  listed  aliovc. 
This  may  also  be  seen  by  inspection  of  Crocco^s 
equations  for  the  laminar  compressible  boundary 
layer  (ref.  22). 


EQUIVALENCE  WITH  PERFECT  GASES 

When  specialized  for  constant-pressure  turbu¬ 
lent  boundary-layer  flows  of  thermally  and 
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calorically  perfect  gases,  the  similarity  parameters 
and  variables  appear  in  the  following  form : 


Rci 


>  Controllable  parameters 


p(f) 

Pr{f) 

1{t) 

m,  v) 


V State  variables 


^Turbulence  variables 


Since  air  is  the  gas  of  principal  interest,  it  is  con¬ 
venient  to  equate  the  Mach  number  parameter, 
(7— for  a  gas  other  than  air,  to  the  param¬ 
eter  in  air  and  to  define  an  equivalent  air  Mach 
number,  Ma. 


Thus,  a  gas,  flowing  at  an  actual  free-stream 
Mach  number,  Mg,  has  an  equivalent  air  Mach 
number,  Ma,  as  defined  by  equation  (9). 

It  is  of  interest  to  consider  the  validity  of  the 
parameter,  {y—l)MJ,  for  flows  other  than 
constant-pressure  turbulent  boundary  layers.  As 
noted  above,  this  parameter  is  valid  for  either 
laminar  or  turbulent  boundary  layers  at  constant 
pressure.  It  is  also  valid  for  comparing  two 
boundary  layers  with  matched  dimensionless  pres¬ 
sure  gradients.  However,  in  this  case  the  dimen 
sionless  differential  equations  must  be  normalized 
with  respect  to  constant  reference  quantities,  not 
to  the  local  free-stream  quantities.  The  param¬ 
eter,  (7— l)Moo^  is  not  valid  for  external  flows. 
This  means  that  in  matching  dimensionless  pres¬ 
sure  gradients  around  two  bodies  in  external 


flows,  the  bodies  must,  in  general,  have  different 
shapes. 

EQUIVALENCE  OF  AIR  AND  HELIUM 

If  air  and  helium  are  considered  perfect  gases  in 
the  constant-pressure  turbulent  boundary  layer, 
the  parameters  and  terms  that  must  be  matched 
for  dynamic  similarity  are  those  listed  in  the 
previous  subsection. 

General  comparison, — The  controllable  param¬ 
eters,  Rei,  {y~~l)MJ,  and  can  easily  be 
matched  between  an  air  and  a  helium  boundary 
layer.  It  is  of  interest  to  note  that  for  adiabatic 
boimdary  layers,  the  boundary  condition  on  Tw 
is  automatically  matched  when  the  Mach  number 
parameter,  (7— 1)M«^  is  matched  for  gases  with 
equal  recovery  factors.  (As  reported  in  appendix 
A,  the  recovery  factors  for  air  and  helium  were 
found  to  be  essentially  the  same.)  The  equiva¬ 
lent  air  Mach  number  for  a  helium  flow  may  be 
written  as  (eq.  (9)): 

(10) 

V  TAlr  ^ 

The  state  variables  that  should  be  matched  are 
the  quantities,  p(T)y  Pr{T),  and  )[i(T).^The 
dimensionless  density,  p(T)  is  equal  to  IjT  for 
both  air  and  helium  when  they  are  considered 
perfect  gases.  The  Prandtl  munber,  Pr{T),  will 
be  approximately  matched  between  air  and  helium 
over  a  considerable  temperature  range  as  shown 
in  figure  1.  The  figure  also  shows  that  the  Prandtl 
number  for  both  air  and  helium  varies  slowly  with 
temperature  over  a  wide  temperature  range  (refs. 
23,  24,  25).  (In  the  regime  of  continuum  flow 
calculations,  the  Prandtl  number  is  often  taken  as 


T,^R 

Figure  1. — Prandtl  numbers  of  air  and  helium. 
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Figure  2.-“Viscosities  of  air  and  helium. 

constant.)  The  dimensionless  viscosit y, 
cannot  be  exactly  matched  between  air  and  lielinm. 
Figure  2  shows  a  logarithmic  plot  of  viscosities 
of  air  and  helium.  The  slopes  of  the  curves  rep¬ 
resent  the  exponents  in  a  power  law  for  viscosity 
and  indicate  the  degree  ot  matching  of  dimension¬ 
less  viscosities. 

For  air  at  high  temperatures: 

For  air  at  low  temperatures: 

The  dimensionless  viscosity  of  helium  over  a  wide 
temperature  range  may  be  closely  approximated  as : 


It  is  of  interest  to  observe  from  the  figure  that  at 
the  higher  Mach  numbers,  tlie  dimensionless  vis¬ 
cosity  of  air  at  ordinary  wind-tunnel  temperatures 
does  not  perfectly  match  air  at  flight  temperatures. 
It  is  also  seen  that  tlie  dimensionless  viscosity  of 
helium  matches  air  at  flight  temperatures  and 
high  Mach  numbers  better  than  does  air  at 
ordinary  wind-tunnel  temperatures. 

As  mentioned  in  the  subsection,  Analysis  of 
Equivalence,  it  is  not  possible  to  know  in  advance 
whether  the  turbulence  terms, 


nrc  actiinlh'  mntcliod  botwocn  an  air  and  a  liclium 
turbulent  boundary  layer.  Ii  dynamic  similarity 
between  an  air  and  a  helium  turbulent  boundary 
la^^er  can  be  demonstrated  experimentally,  it 
should  be  concluded  that  the  turbulence  terms  arc 
sufhciently  well  matched. 

Specific  comparison. — In  attempting  to  set  up 
equivalent  air  and  helium  boundary  layers  in  a 
wind  tunnel,  at  ordinary  wind-tunnel  temperature 
(e.g.,  R),  one  quantity  that  will  have  an 

appreciable  mismatch  in  a  portion  of  the  boundary 
layer  is  the  dimensionless  viscosity,  This 

mismatch  is  most  serious  at  the  wall  (e.g.,  with 
and  ^^==000°  R, 1.5  When 

the  dimensionless  equations  of  motion  (4),  (5), 
and  (6)  are  considered,  it  is  seen  tliat  a  mismatch 
in  Jl(T)  must  have  some  effect  on  the  other  quan¬ 
tities,  but  the  magnitude  of  the  effect  is  difficult 
to  estimate.  (In  ref.  20  it  is  predicted  that  a 
mismatch  in  Ji  will  have  a  small  effect  on  turbulent 
skin  friction,  but  a  large  effect  on  laminar  skin 
friction.)  In  equations  (5)  and  (6)  the  product 
(l//?cdM  occurs.  This  can  be  written 


{llRer^  =  {^l^,ypAyne„)  (11) 


where  7?c„. 


is  the  “vv^all  property”  Reyn¬ 


olds  number.  The  product  (l/7iV/)pL  can  bo 
matched  in  any  selected  portion  of  the  boundary 
layer  b}^  mismatching  Rci  to  compensat{^  for  any 
mismatch  in  Ju.  It  has  been  thought  that  (I/Ac^m 
should  be  matched  somewlu're  in  the  sublayer 
(possibly  at  the  wall)  as  this  sliould  be  the  region 
wliere  viscosity  is  most  important.  (Away  horn 
the  wall  €  becomes  much  larger  than  Ji  so  that-  ju 
should  be  less  important  away  fjom  the  wall.) 

The  method  of  matching  (1/Jiei)fl  near  or  at  the 
wall  is  easily  seen  by  referring  to  equation  (11). 
Near  the  wall  /i/mu-  will  be  nearly  unity,  hence 
p/jjijr  will  be  reasonably  well  matched  regardl(‘ss  ol 
the  gas.  The  dimensionless  density  p„;ismatclied 
between  two  flows  being  compared  at  the  same 
and  It  is,  then,  necessary  to  match 

in  order  to  secure  a  good  match  of  (l/7fcdM 
or  at  the  wall.  This  results  in  a  poor  match  of 
(l/7iVdju  away  from  the  wall  (where  p  is  possibly 
not  important).  The  practical  question  of  whether 
{llEei)jl  should  better  be  matched  near  the  wall 
or  in  the  free  stream  can  be  best  answered  ex¬ 
perimentally  by  determining  whether  measured 
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Cf  values  correlate  with  Bei  or  Bcy,  (it  either). 
(Obviously,  if  there  were  no  mismatch  in  Jx,  this 
question  would  not  arise  and  it  would  then  be 
immaterial  whether  Rei  or  Ecy,  were  used  in  the 
comparisons.) 

In  determining  the  practical  equivalence  of 
turbulent  boundary  layers  in  air  and  helium,  it 
will  be  useful  to  consider  boundary-layer  profiles 
with  the  variables  in  several  forms.  The  dimen¬ 
sionless  velocity,  S,  may  be  written  as : 

^  M  1'^ 

This  means  that  M/Mcx,  will  be  equal  at  corre¬ 
sponding  points  in  two  perfectly  equivalent 
boundary  layers.  (The  M  values  themselves  in 
the  air  and  helium  boundary  layers  will  be  in  the 
ratio  —1.29  at  corresponding  points.)  Also,  cor¬ 
responding  points  in  two  perfectly  matched  bound¬ 
ary  layers  are  at  equal  values  of  y/O  which  can  be 
shown  as  follows: 


y  vl 

7)1 

e  0 

Cf{xi)dxi 

'o 

y 

V 

e 

(1/2) 

'o 

The  relations  above  are  for  perfectly  matched 
boundary  layers.  In  view  of  the  mismatch  in  Ji 
between  air  and  helium  at  wind-tunnel  tempera¬ 
tures,  it  is  desirable  to  compare  experimental 
profiles  to  assess  tlie  effect  of  the  Ju  mismatch  and 
to  assist  in  determining  whether  equivalence  can 
be  considered  to  exist  for  practical  purposes. 

Profile  comparisons  can  also  be  made  using 
vv 

'ii^—ujur  and  as  coordinates,  with  the 

friction  velocity,  Ur,  defined  as: 


The  quantities,  u*  and  as  defined,  are  nor¬ 
malized  with  respect  to  wall  properties.  The 
equations  for  and  ?/*  can  be  written  as: 


±  1^  /I 
'A  2Vf. 


From  this  it  is  clear  that  if  two  boundary  layers 
are  perfectly  equivalent  (having  equal  values  of 
Gf  and  Ty,  and  equal  values  of  u  at  corresponding 
points),  the  quantity,  7^*,  will  be  matched  at  corre¬ 
sponding  points.  The  quantity  y*  will  be  matched 
at  corresponding  points  if  Cf  is  shown  to  correlate 
with  Rcw  for  two  boundary  layers.  On  the  other 
hand,  if  Cf  correlates  with  Rei  for  two  boundary 
layers,  then  the  i/*  values  between  air  and  helium 
boundary  layers  will  not  be  matched  at  corre¬ 
sponding  points  (at  wind-tunnel  temperatures) 
because  of  the  mismatch  in  iiy,.  At  a  total  tem¬ 
perature  of  600°  R  and  an  Ma  value  of  4.2  (giving 
mismatch  in  y*  values  will 

be  y*He~1.5  ?/*Air  for  corresponding  points.  Ex¬ 
perimental  profiles  showing  the  variables  in  the 
several  forms  mentioned  above  should  (in  addition 
to  skin-friction  measurements)  assist  in  deter¬ 
mining  whether  it  is  better  to  compare  on  an  Rei 
or  an  Rew  basis. 

ADVANTAGES  OF  USING  HELIUM 

As  IS  well  known,  it  is  difficult  to  obtain  high 
Mach  numbers  and  high  Reynolds  numbers  simul¬ 
taneously  with  air  in  a  wind  tunnel.  In  order  to 
avoid  condensation,  it  is  necessary  to  heat  the  air, 
with  a  resulting  loss  in  density  and  Reynolds 
number.  It  seems  desirable,  then,  to  use  a  gas 
such  as  helium  that  does  not  condense  at  the  high 
Mach  numbers  of  interest.  As  pointed  out  earlier, 
helium  under  wind-tunnel  conditions  provides  a 
good  simulation  of  air  under  flight  conditions  pro¬ 
vided  imperfect  gas  effects  are  not  significant. 
As  shown  in  figure  3,  for  a  given  total  pressure, 
higher  Reynolds  numbers  can  be  obtained  at 
hypersonic  Mach  numbers  with  helium  than  with 
air  in  a  wind  tunnel.  As  an  example,  at  an  equiva¬ 
lent  air  Mach  number  of  10,  the  Reynolds  number 
factor  in  favor  ol  helium  is  about  13  to  1.  (This 
factor  has  a  slow  dependence  on  pressure  when 
the  total  temperature  of  air  is  taken  as  the  mini¬ 
mum  necessary  to  avoid  condensation  at  a  given 
Mach  number;  in  figure  3  a  total  pressure  of  1000 
psi  was  used.)  This  advantage  of  helium  is  due 
to  several  factors,  the  principal  one  being  that  no 
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Ficjure  3. — Variation  of  Reynolds  niinilx'r  }>or  foot  with 
cqnivah'nt  air  Mach  number  for  air  and  lielium; 

1000  psi. 


heating  of  helium  is  required.  Otlier  factors 
that  contribute  to  tliis  result  are  the  different 
ratios  of  specific  heats,  and  the  fact  that  with 
equal  equivalent  air  Mach  numbers,  testing  in 
helium  can  be  conducted  at  a  lower  actual  Alach 
number. 

APPARATUS  AND  TEST  CONDITIONS 

THE  1-  BY  10-INCII  BOUNDARY-LAYER  CHANNEF. 

The  1-  by  10-inch  boundary-layer  channel  was 
designed  to  provide  high  Mach  numbers  and  high 
Reynolds  numbers  simultaneousiy.  In  figure  4 
the  test  domain  of  the  boundary-layer  channel  is 
illustrated  in  terms  of  Mach  number  and  Reymolds 


Figure  4. — Test  domain  of  the  1-  by  10-inch  boundary 
layer  channel. 


number.  For  purposes  of  comparison,  also  shown 
is  the  so-called  corridor  of  continuous  flight,  which 
is  the  approximate  Mach  number  Reynolds  num¬ 
ber  range  lor  continuous,  sustained  flight  of  air¬ 
planelike  vehicles.  (A  characteristic  length  of  50 
feet  was  used  in  the  calculation;  the  upper  bound 
of  the  corridor  was  taken  at  an  equilibrium  sur¬ 
face  temperature  of  2000°  F,  and  the  lower  bound 
was  estimated  using  a  dynamic  pressure  of  100 
psf.)  Most  previous  skin-friction  data  were  taken 
in  the  area  to  the  left  and  below  the  lower  bound¬ 
ary  of  the  corridor  of  continuous  flight.  It  is 
seen  tliat  the  boundary-layer  channel  has  extended 
the  Mach  number  Reynolds  number  range  into  a 
regime  of  considerable  interest.  In  order  to  cover 
the  range  shown,  air  was  used  as  the  working  fluid 
up  to  a  Alach  number  of  4.2,  and  helium  was 
used  from  equivalent  air  Mach  numbers  of  4.2 
up  to  9.9. 

A  sketch  of  the  equipment  layout  is  shown  in 
figure  5;  figure  6  is  a  schematic  drawing  of  the 
test  region;  and  figure  7  is  a  photograph  of  tlie 
equipment.  Altliough  this  apparatus  is  called  a 
channel  (because  of  its  dimensions),  the  series  of 
tests  being  reported  were  not  for  channel  flows. 
During  the  present  series  of  tests,  there  was  always 
a  core  of  potential  flow  between  the  top  and  the 
bottom  boundary  layers. 

The  1-  by  10 -in eh  boundary-layer  channel  is  a 
blowdown  type  wind  tunnel.  The  gas  used  as  a 
working  fluid  is  stored  at  2400  psi  in  two  40,000 
cubic  foot  capacity  gas  trailers.  As  shown  in 
figure  5  tlie  gas  flows  from  the  trailers  through 
the  inlet  line  into  the  heat  exchanger,  and  then 
into  the  settling  chamber  of  the  tunnel  at  the 
desired  temperature  and  ju’essure.  The  total  pres¬ 
sure  (measured  in  the  settling  chamber)  was  con¬ 
trolled  by  valves  near  the  trailers.  The  total 
temperature  was  controlled  by  means  of  a  heat 
exchanger  consisting  of  a  thermal  mass  of  3,000 
pounds  of  copper  tubes.  All  heating  of  the  tlier- 
mal  mass  was  performed  prior  to  running  in  order 
to  preset  the  desired  total  temperatures.  Runs 
were  terminated  wlien  the  heat  exchanger  could 
no  longer  maintain  the  preset  total  temperature. 
In  the  present  series  of  tests  the  total  temperat.ure 
was  preset  relative  to  the  ambient  wall  tempera¬ 
ture  in  the  tunnel  test  section  so  that  the  boundary 
layers  were  adiabatic  at  the  walls.  Tliis  procedure 
gave  total  temperatures  for  all  runs  of  approxi¬ 
mately  600°  R.  (Sec  appendix  A.) 
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High-pressure  gos  from  trailers 


Figure  6. — Details  of  test  region  of  boundary-layer 
channel. 


The  tunnel  dimensions  (fig.  6)  were  designed  to 
give  an  essentially  two-dimensional  flow.  The 
top  plate  of  the  tunnel  is  completely  flat;  the 
nozzle  contouring  is  entirely  on  the  bottom  wall. 
Four  interchangeable  nozzles  were  used,  three 
contoured  for  supersonic  flows  and  one  (a  flat 
plate)  for  subsonic  flows.  Mach  number  and 
pressure  gradient  settings  were  made  by  the  selec¬ 
tion  and  adjustment  of  the  interchangeable  nozzle 
blocks.  Keynolds  number  control  was  obtained 


by  varying  the  total  pressure  for  a  given  tunnel 
geometry. 

As  shown  in  figure  5,  the  apparatus  has  two 
exhaust  routes  available,  atmospheric  and  vacuum. 
During  supersonic  running,  the  exhaust  was  al¬ 
ways  to  the  vacuum  exhaust  system.  An  opti¬ 
mum  compression  ratio  was  achieved  by  means 
of  an  adjustable  second  thioat.  The  atmospheric 
exhaust  was  used  for  subsonic  running.  For  the 
subsonic  runs,  Mach  number  control  was  obtained 
with  the  (supersonic)  diffuser  throat  operating  as 
a  sonic  choke  downstream  of  the  test  section. 

In  the  series  of  tests  conducted,  the  Mach  num¬ 
ber  and  Reynolds  number  were  varied  to  the  limits 
of  the  present  capability  of  the  equipment.  The 
upper  limit  of  the  Mach  number  (and  the  lower 
limit  of  the  Reynolds  number  at  a  given  Mach 
number)  were  determined  by  the  exhaust  system 
back  pressure  (compression  ratio  available) ,  while 
the  maximum  Reynolds  numbers  were  dictated  by 
safety  considerations.  The  maximum  settling 
chamber  pressure  was  700  psig  for  supersonic  flows 
and  310  psig  for  subsonic  flows. 
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Figure  7. — The  1-  by  lO-inch  boundary-layer  channel  witliout  instruinenlation. 


SKIN-FRICTION  ELEMENT 

Local  skill  friction  was  measured  directly  by 
means  of  the  skin-friction  element  shown  in  figures 
8  and  9.  The  element  consisted  of  u  disk  2}i 
inches  in  diameter  suspended  from  two  flexures. 
The  sti’oamwise  force  on  the  disk  was  mcuisured  by 
means  of  a  difl’erential  transformer.  The  trans¬ 
former  itself  maintained  a  fixed  position  while  two 
iron  magnetic  return  pieces  on  the  disk  completed 
the  magnetic  circuit.  All  of  the  metal  parts  of  the 
element  not  in  the  magnetic  circuit  wei’c  made  of 
K-monel.  In  the  neutral  position,  the  gap  be¬ 
tween  the  disk  and  its  housing  was  0.005  inch. 
The  element  had  no  nulling  device  and  simply 
deflected  under  load. 
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The  top  face  of  the  disk  was  formed  as  a  dash 
pot  which  contained  1000  centistoke  silicone  oil  as 
a  damping  fluid.  The  disk  in  its  mounting  was 
electrically  insulated  from  the  element  housing  so 
that  an  electric  fouling  meter  could  be  used  to 
indicate  contact  between  the  disk  and  the  housins:. 
Tlie  mounts  for  the  flexures  constrained  the  ends 
of  the  flexures  to  remain  pai’allel  and  vertical. 
This  meant  that  the  disk  was  constrained  to  re¬ 
main  parallel  to  its  no-load  position  (horizontal). 
Interchangeable  pairs  of  flexures  were  calibrated 
and  used  for  different  loading  ranges.  The 
approximate  force  at  a  full  load  deflection  of  ap- 
proximatel}^  0.005  inch  for  tlie  pair  of  flexures  used 
for  most  of  the  tests  was  0.1  pound.  The  disk 
was  determined  to  be  flush  with  its  housing  to 
within  10  microinches  by  means  of  an  optical 
interferometer.  Subsecpient  tests  to  determine 
fliishness  requirements  showed  that  for  a  Mach 
number  of  4.2,  the  element  disk  could  be  displaced 
as  much  as  200  microinches  in  or  out  of  the  housimr 
with  no  noticeable  change  in  skin  friction.  The 
surface  of  the  disk  was  finished  to  2  microinches, 
the  housing  to  4  microinches,  while  the  top  plate 
of  the  channel  was  finished  to  15  microinches. 
These  finishes  were  selected  to  provide  “hydraulic¬ 
ally  smooth’^  surfaces  over  the  range  of  variables 
tested.  The  finishes  on  the  disk  and  housing  were 


measured  by  means  of  optical  flats  and  on  the  top 
plate  by  means  of  a  profilometer.  Twelve  static 
orifices  in  the  gap  surface  of  the  housing  surround¬ 
ing  the  disk  provided  pressure  data  for  calculating 
tl)e  buoyancy  force  on  the  disk.  These  pressures 
were  read  on  a  small  manometer  boairl  neai*  the 
element  which  used  10  centistoke  silicone  oil. 

The  output  of  the  difTerential  transformer  was 
fed  into  a  bridge  circuit  sei  vo  witli  a  digital  i-ead- 
oiit.  The  electrical  circuitry  is  shown  in  figure  10. 
The  calibration  of  the  element  and  servo  combina¬ 
tion  showed  a  repeatability  of  0.1  percent  of  the 
full  load  reading  and  was  nearly  linear  over  tlie 
full  range  (depai'ted  approximately  2  percent  fi’om 
linearity  at  full  load).  Over  the  range  of  tem¬ 
peratures  used,  the  temperature  effect  on  the 
calibration  was  almost  negligilile.  (Calibration 
curves  taken  at  40°  F  and  83°  F  showed  a  slope 
difference  of  approximately  %  percent.)  A  descrip¬ 
tion  of  the  electrical  circuitry  of  the  element  and 
servo  is  given  in  appendix  B. 

GENERAL  INSTRUMENTATION 

Static-pressure  distributions  were  obtained 
through  orifices  in  the  top  plate  of  the  tunnel  (fig. 
11).  The  orifices  were  connected  to  a  30-tubc 
manometer  board  that  used  10  centistoke  silicone 
oil  and  was  referenced  to  a  near  vacuum.  The 


Figure  10. — Circuitry  of  skin-friction  element  and  servo. 
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total  pressure  was  measured  at  the  settling  cham¬ 
ber  with  a  Bourdon  gage.  Impact  pressures  in 
the  potential  flow  core  were  measmed  by  two  fixed 
probes  connected  to  pressure  transducers  which 
were  electrically  connected  to  bridge-type  servos 
with  digital  readouts.  Boundary-layer  impact 
pressures  were  measured  by  two  movable  probes 
connected  to  pressure  transducers,  and  located 
1%  inches  on  either  side  of  the  element  center  line. 
The  boundary-layer  probes  were  flattened  tubes 
with  an  inside  height  of  0.002  inch,  an  outside 
height  of  0.004  inch,  and  a  width  of  0.045  inch. 
The  traversing  of  the  probes  was  controlled  re¬ 
motely  and  could  be  positioned  to  within  0.002 
inch. 

Temperatures  were  measured  at  17  locations 
including  the  settling  chamber,  the  skin-friction 
element,  several  locations  and  depths  in  the  top 
plate,  the  manometer  boards,  and  the  heat  ex¬ 
changer.  Iron  constant  an  thermocouples  with 
self-balancing  potentiometers  were  used.  The 
settling  chamber  temperature  was  continuously 
recorded. 

PROCEDURE  AND  DATA  REDUCTION 

Tests  conducted  were  of  three  types:  skin- 
friction  tests,  transition  tests  to  determine  the 
virtual  origin  of  the  turbulent  boundary  layer, 
and  boundary-layer  surveys. 


SKIN-FRlCTION  TESTS  (AIR  AND  HELIUM) 

Skin-friction  tests  in  air  were  performed  at  both 
subsonic  and  supersonic  Mach  numbers  while  the 
tests  using  helium  as  the  working  fluid  were  made 
at  supersonic  Mach  numbers  only.  The  quan¬ 
tities  directly  measured  were :  the  force  on  the  skin- 
friction  element,  the  total  pressure  (in  the  settling 
chamber),  impact  pressures  (usually  2)  in  the 
potential  flow  core  near  the  element,  the  static 
pressure  at  the  skin-friction  element,  12  buoyancy 
pressures  around  the  element  disk,  30  static  pres¬ 
sures  on  the  top  plate  of  the  test  section,  the  total 
temperature  (in  the  settling  chamber),  and  other 
temperatures  of  interest. 

The  free-stream  Mach  number,  M,»,  and  the 
dynamic  pressure,  g«>,  were  obtained  from  the 
measured  static  and  total  pressure  and  the  usual 
isentropic  flow  equations.  The  exisrence  of  an 
isentropic  core  was  verified  by  a  comparison  of 
impact  pressures  calculated  from  isentropic  and 
shock  relations  with  measured  impact  pressures. 
These  pressures  agreed  within  4  percent. 

In  the  determination  of  Eeynolds  number,  the 
length  used  was  from  the  virtual  origin  of  the 
turbulent  boundary  layer.  The  determination  of 
the  virtual  origin  is  discussed  in  the  following 
subsection.  For  the  viscosity  of  air  (fig.  2), 
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Sutherland’s  formula  in  tlie  following  form  was 
used  (refs.  27  and  28): 

7^3/2 

2.270  soc/ft^’  (12) 

For  helium,  in  the  temperature  range  at  whi(‘li 
tests  were  run,  the  viscosity-temperature  relation¬ 
ship  used  was  (hg.  2) : 

M-7.35T«-«-'X10-®  lb  sec/ft^  (13) 

This  formula  is  based  on  a  faired  curve  through 
available  experimental  data  (ref.  24). 

The  buoyancy  force  on  the  disk  of  the  skin- 
friction  element  was  obtained  by  numei*ically 
integrating  the  sti’eamwise  component  of  forces 
due  to  the  pressures  acting  on  the  eflective 
projected  area  of  the  disk.  The  following  sketch 
shows  the  disk  dimensions  and  orifice  locations: 


Ficuhe  12. — IP'lnlivc  map:iiitud(‘  of  buoyniu'v  forco  on 
skiii-frici  ion  element. 


Tests  at  a  given  Mach  number  were  conducted 
witli  a  fixed  nozzle  and  test-section  geometry. 
Consequently,  there  was  a  slight  variation  in 
Mach  number  with  Reynolds  number  as  a  result 
of  tlie  variation  in  [)Oundary-layer  thickness  as 
shown  in  figure  13.  In  order  to  present  skin- 


Defining  the  buo^^ancy  force,  F^uovj  as  positive 
in  the  upstream  direction,  and  defining  Ftomi 
as  the  total  force  sensed  by  the  disk  of  the  skin- 
friction  element,  one  can  write: 


W^  —  1^  Skin  friction  totnl~\~  Fi 


buoy 


Reynolds  number,  millions 


Fioxtre  1:L — Variation  of  txmiH'l  Mach  mimbnr  willi  Rxyn- 
olds  inimb('r  for  various  nozzle  settin^r.s. 


where  A  is  the  flat  surface  area  of  the  disk. 
The  experimentally  obtained  ratio  of  the  buoy¬ 
ancy  force  to  the  skin-friction  force,  FtuoylFskin 
friction,  as  a  function  of  Reynolds  number  for  the 
various  Mach  numbers  used  is  shown  plotted  in 
figure  12.  It  is  seen  that  the  buoyancy  forces 
were  relatively  small  except  for  the  case  with 
iV/a=6.7  (see  subsection  below  entitled  Accuracy). 
The  coefficient  of  local  skin  friction,  Cf^  is  given 
by  the  defining  equation: 


friction  results  at  fixed  values  of  Mach  mimber, 
the  measured  skin-friction  coeffidents  were  cor¬ 
rected  by  means  of  tlie  following  ex])ression: 


Althougli  the  factor  (A^V/AMco)  could  have  been 
determined  by  experiment,  it  was  more  conven¬ 
iently  determined  by  the  T'  method  (ref.  19). 
Any  errors  introduced  by  this  procedure  arc 
unimportant  since  the  corrections  were  always 
less  than  1  percent. 
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TRANSITION  TESTS  TO  DETERMINE  THE  VIRTUAL  ORIGIN 

As  noted  previously,  the  length  used  in  the 
calculation  of  Reynolds  numbers  was  taken  from 
the  approximate  location  of  the  virtual  origin  of 
the  turbulent  boundary  layer  to  the  center  of 
the  skin-friction  element.  In  locating  the  virtual 
origin  for  supersonic  flows,  it  was  elected  to 
determine  the  point  of  maximum  local  skin  fric¬ 
tion,  which  is  known  to  be  near  the  point  of 
transition  (ref.  29).  The  point  of  maximum 
local  skin  friction  was  determined  by  means  of 
impact  probes  at  the  wall.  These  probes  were 
hollow  needles  of  0.020  inch  diameter  with  a 
slot  of  height  0,002  inch  for  pressure-sensing. 
The  needles  were  fitted  in  the  top  plate  static 
orifices  and  were  adjusted  so  that  the  slot  just 
cleared  the  wall  surface.  The  impact  probes  at 
the  wall  were  located  in  fixed  positions  on  the  top 
plate  and  the  total  pressure  was  varied.  It  can 
be  shown  that: 

where  Ap  is  the  difference  between  the  impact 
pressure  sensed  by  the  probe  and  the  local  static 
pressure,  pt  is  the  settling  chamber  pressure,  and 
k  is  a  proportionality  constant.  This  crude 
approximation  was  adequate  for  the  purpose  used ; 
it  was  not  necessary  to  calibrate  the  probes,  since 
the  only  interest  was  in  knowing  whether  the 
boundary  layer  was  laminar,  turbulent,  or  in 
transition.  For  each  x-station  at  which  the  probes 
were  placed,  a  plot  was  obtained  as  illustrated  in 
the  following  sketch: 


From  the  group  of  plots  typified  by  the  sketch 
above,  a  plot  was  made  (for  each  Mach  number) 
showing  the  location  of  the  point  of  maximum 
local  skin-friction  as  a  function  of  the  total 
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pressure.  This  type  of  plot  is  illustrated  in  the 
following  sketch: 


With  decreasing  total  pressure,  pt,  the  point  of 
transition  moves  downstream  and  eventually 
reaches  the  element  itself.  Data  were  not  reduced 
or  reported  where  transition  was  completed  at  a 
distance  closer  than  8  inches  upstream  of  the 
element’s  center.  This  distance  corresponded  to 
approximately  4  disk  diameters  or  to  a  minimum 
of  approximately  12  boundary-layer  thicknesses. 
This  distance  limitation  was  selected  because  it 
was  felt  that  the  turbulent  boundary  layer  was 
probably  not  in  equilibrium  (free  from  historical 
effects)  over  shorter  distances.  Another  factor 
in  selecting  an  8  inch  minimum  length  of  run  was 
that  the  disk  diameter  was  required  to  be  small 
relative  to  the  length  of  run  of  the  turbulent 
boundary  layer,  inasmuch  as  the  coefficient  of 
local  skin  friction  as  calculated  and  reported  waS 
actually  an  average  value  over  the  disk  area. 

In  order  to  have  a  convenient  method  for  con¬ 
verting  from  Rbl  to  Be^,  the  data  typified  by  the 
preceding  sketch  were  replotted  into  a  working 
plot  in  figure  14  with  Bcl  as  abscissa,  and  BexjBeL 
as  ordinate.  The  ordinate  is  thus  a  factor  to  be 
multiplied  by  Bei,^  to  obtain  the  Reynolds  number 
measured  from  the  virtual  origin,  Be^.  The 
dotted  extensions  of  the  curves  are  extrapolations 
based  on  a  constant  Reynolds  number  of  transi¬ 
tion  which  was  taken  as  the  value  obtained  at  the 
farthest  upstream  station  at  which  measurements 
were  made.  These  Reynolds  numbers  are  as 
indicated  in  the  figure.  It  is  seen  in  figure  14 
that  the  Reynolds  number  factor  is  nearly  unity 
at  the  higher  Reynolds  numbers.  This  means 
that  with  supersonic  flow  at  high  Reynolds 
numbers  the  virtual  origin  was  close  to  the  nozzle 
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Re^  ,  minions 


Figure  14. — Factor  for  determining  Heynolds  number 
measured  from  the  virtual  origin. 

throat.  Except  at  il/a=9.9,  the  great  bulk  of 
skin-friction  data  were  obtained  at  Reynolds 
numbers  greater  than  shown  on  figure  14.  For 
these  cases  Bcx  was  obtained  by  subtracting  from 
Bcl  the  Reynolds  number  of  transition  (as  labeled 
at  the  tops  of  the  curves). 

With  subsonic  fiow^  a  different  procedure  was 
used  for  obtaining  the  virtual  origin.  The  probe 
technique  described  above  showed  that,  at  the 
Reynolds  numbers  of  interest,  transition  occurred 
ahead  of  the  orifice  farthest  upstream  in  the  top 
plate.  For  calculations  of  Reynolds  number,  the 
point  of  transition  was  taken  to  be  the  junction 
of  the  sharply  converging  transition  section  (fig. 
6)  and  the  flat  plate  section.  This  point  was  2.60 
inches  upstream  of  what  would  have  been  the 
location  of  the  nozzle  throat  for  supersonic  runs 
(fig.  11). 

BOUNDARY-LAYER  SURVEYS 

Boundary -layer  surveys  at  the  longitudinal 
position  of  the  skin-friction  element  were  obtained 


at  eacli  of  the  supersonic  Mach  numbers  and  at 
selected  Reynolds  numbers.  The  surveys  were 
performed  by  the  usual  method  of  ])itot  travei*sing. 
Fi’cqiiently  two  sui‘veys  were  taken  simulta¬ 
neously,  one  on  each  side  of  the  skin-friction 
element.  The  momeTitum  thickness,  0,  was  of)- 
tained  from  the  defining  equation  (see  symbols) 
by  graphical  integration.  Constant  static  pres- 
sui'e,  constant  total  temperature  (settling  cham})er 
temperature),  and  constant  specific  heats  across 
the  boundary  layer  were  assumed  in  the  calcula¬ 
tion  of  u  and  p  through  the  boundary  layer  and 

of  6,  In  the  calculations  of  and 

the  friction  velocity  was  calculated  as  Ur~-\ 

V  Pto 

where  was  measured  by  the  skin-friction 
element. 

ACCURACY 

All  effort  was  made  to  use  procedures  tliat  would 
give  the  best  possible  accuracy  of  data  from  the 
equipment.  To  keep  random  errors  witliin  small 
bounds,  skin-friction  measurements  with  air  were 
repeated  five  times,  and  data  with  helium,  three 
times.  The  scatter  of  the  Cf  data  indicated  a 
probable  random  error  of  approximately  2  percent. 
Possible  sources  of  systematic  errors  arose  from 
the  correction  for  the  buoyancy  force,  from  the 
effects  of  a  small  or  irregular  pressure  gradient, 
and  in  the  calculation  of  Reynolds  number.  The 
contributions  of  these  systematic  errors  will  be 
considered  in  the  estimate  of  the  over-all  accuracy 
oiCf. 

An  error  in  the  buoyancy  force  correction  could 
have  been  introduced  through  not  knowing 
exactly  the  effective  area  on  which  the  buoyancy 
pressures  acted,  and  in  not  knowing  the  buoyancy 
pressure  over  every  infinitesimal  element  of  area. 
It  is  estimated  that  the  buoyancy  force  is  in  error 
by  less  than  10  percent.  Since  the  buoyancy 
force  was  generally  less  than  5  percent  of  the  skin- 
friction  force,  and  never  more  than  22  percent 
(fig.  12),  tlm  error  in  the  skin-friction  coefficient 
due  to  buoyancy  force  error  is  probably  less  than 
2  percent.  As  a  check  on  this,  a  number  of  skin- 
friction  tests  were  run  at  a  Mach  number  of  4.2 
using  a  disk  in  the  skin-friction  element  with  a  side 
projected  area  of  four  times  that  used  in  the  regiilar 
tests.  This  disk  gave  buoyancy  forces  that  were 
approximately  20  percent  of  the  skin-friction 
force,  but  the  skin-friction  coefficients  obtained 
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with  the  two  disks  checked  each  other  within  2 
percent.  From  this  experience,  it  is  believed  that 
the  buoyancy  forces  at  a  Mach  number  of  6.7  were 
not  excessive.  (The  irregularity  of  the  buoyancy 
force  trend  at  this  Mach  number  shown  in  figure 
12  is  believed  to  be  due  to  a  slight  wave  pattern 
in  the  flow.) 

The  deviations  of  static  pressure  from  zero 
pressure  gradient  as  obtained  in  typical  tests  are 
shown  in  figure  15.  The  small  wave  pattern  at 
it/a = 6.7  can  be  seen.  Also  shown  is  a  small  favor¬ 
able  pressure  gradient  at  it/a—2.95.  Deviations 


Air  /l^^w2.95 
Helium  4.2 

Air  4.2 

Helium  6,7 

Helium  9.9 


20x106) 


(a)  Longitudinal  pressure  distributions. 


Distonce  from  tunnel  ^  ,  inches' 


(b)  Transverse  pressure  distributions  at  various  distances, 
L,  from  nozzle  throat. 

Figure  15. — Typical  static-pressure  distributions  for 

26X108. 


from  zero  pressure  gradient  of  these  amounts  were 
used  in  a  special  series  ol  tests  at  Ma— 4.2  in  order 
to  determine  the  eflPect  of  small  pressure  gradients. 
It  was  found  that  the  change  in  Cf  from  the  zero 
pressure  gradient  case  was  less  than  the  small 
scatter  of  the  data.  It  is  believed  that  aU  of  the 
data  have  a  probable  error  in  Of  due  to  pressure 
gradient  of  less  than  2  percent. 

Errors  in  the  calculation  of  Reynolds  number 
would  have  been  almost  entirely  due  to  errors  in 
estimating  the  location  of  the  virtual  origin.  The 
true  location  of  the  virtual  origin  was  probably 
slightly  upstream  of  the  point  of  maximum  Gfy 
which  was  used  for  supersonic  flows.  It  seems 
almost  certain  that  the  virtual  origin  for  super¬ 
sonic  flows  was  never  upstream  of  the  nozzle 
throat.  (In  any  case  the  boundary  layer  had  to 
be  very  thin  at  the  nozzle  throat.)  With  super¬ 
sonic  flows  at  the  highest  Reynolds  numbers,  the 
point  of  maximum  Of  (except  at  was 

close  to  the  nozzle  throat  (fig.  14).  It  is  believed 
that  the  probable  error  in  Reynolds  number  must 
have  been  appreciably  less  than  the  quantity, 
l  —  iReJReL)  (fig.  14).  The  probable  error  in 
Reynolds  number,  then,  is  small  for  data  at 
Reynolds  numbers  near  the  upper  ends  of  the 
curves  in  figure  14  (except  at  A/a=9.9).  As  noted 
previously  the  bulk  of  skin-friction  data  were 
obtained  at  Reynolds  numbers  greater  than  shown 
in  figure  14  (except  at  A/a=9.9),  so  that  probable 
error  in  Reynolds  number  for  all  these  data  must 
have  been  small.  With  decreasing  Reynolds  num¬ 
bers  the  quantity,  l  —  (ReJReL)y  and  hence  the 
probable  error  becomes  progressively  greater. 
Data  with  a  turbulent  boundary-layer  length 
between  8  and  14  inches  (0.34<i?e^//?ei<0.60) 
are  shown  with  flagged  symbols  in  the  skin-friction 
plots  to  denote  that  these  data  are  not  as  accurate 
as  the  higher  Reynolds  number  data.  As  men¬ 
tioned  previously,  data  with  a  turbulent  boundary 
length  of  less  than  8  inches  (i?ea;//?6i<0.34)  are  not 
reported.  At  Mach  numbers  other  than  9.9  it  is 
estimated  that  for  the  unflagged  data  (ReJRezO> 
0.60)  the  probable  errors  in  Reynolds  number  are 
less  than  6  percent  at  the  highest  Reynolds  num¬ 
bers  and  less  than  25  percent  at  the  lowest  Rey¬ 
nolds  numbers.  At  Ma=9.9,  these  probable 
errors  are  estimated  to  be  less  than  15  and  25 
percent,  respectively.  Since  the  slope  of  a  plot 
of  log  Of  against  log  Re^  is  approximately  —1/5, 
these  errors  correspond  approximately  to  errors  in 
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CfOl  1  and  5  percent  at  Macii  numbers  other  than 
9,9,  and  3  and  5  percent  at  M„==9.9. 

Witli  subsonic  flows,  at  the  Reynolds  numbers 
tested,  the  virtual  origin  was  upstream  of  the 
farthest  upstream  orifice  in  the  top  plate.  It  is 
believed  tiiat  the  virtual  origin  was  not  upstream 
of  the  point  assumed  (beginning  of  the  flat-plate 
section)  by  any  large  amount  since  this  would 
have  put  the  virtual  origin  in  a  rapidly  converging 
transition  section  (fig.  6).  It  is  estimated  that 
the  probable  error  in  Reynolds  number  for  sub¬ 
sonic  flows  was  less  than  15  percent.  This  would 
cori-espond  to  a  probable  error  in  Cf  of  approxi¬ 
mately  3  percent. 

ft  is  estimated  that,  for  all  Mach  numbers,  the 
total  probable  error  in  the  measured  skin-friction 
coefficients  is  less  than  5  percent  at  the  highest 
Reynolds  numbers  and  less  than  6  percent  at  the 
lowest  unflagged  Reynolds  numbers.  No  attempt 
has  been  made  to  estimate  the  accuracy  of  the 
flagged  points.  With  decreasing  Reynolds  num¬ 
bers,  these  points  become  progressively  less  accu¬ 
rate  and  tend  to  become  erratic,  as  noted  below  in 
the  section,  TEST  RESULTS  AND  DISCUS¬ 
SION.  The  fact  that  tlio  best  accuracy  is  obtained 
at  the  highest  Reynolds  numbers  is  probably 
fortunate,  since  the  highest  Reynolds  numbers  are 
in  the  regime  of  greatest  interest.. 

Probable  errors  in  boundary-layer  thickness  or 
momentum  thickness  are  estimated  to  be  equiva¬ 
lent  to  errors  in  Re.ynolds  number  and  thus  to  have 
approximately  the  same  probable  errors  as  Reyn¬ 
olds  number.  The  probable  error  in  6,  then, 
becomes  greater  at  the  lower  Reynolds  numbers. 
In  some  cases  simultaneous  profiles  were  talam 
on  either  side  of  the  skin-friction  element.  Span- 
wise  variations  in  6  averaged  approximate! v  12 
percent  and  it  is  believed  that  transverse  varia¬ 
tions  in  the  virtual  origin  location  existed  across 
the  test  section  width.  The  maximum  s])anwise 
variation  in  6  was  14  percent  and  occurred  at  a 
low  Reynolds  number  with  a  I'apidly  shifting 
virtual  origin  (fig.  14). 

TEST  RESULTS  AND  DISCUSSION 

The  data  presented  show  the  effects  of  Mach 
number  and  Reynolds  number  on  local  skin 
friction  with  adiabatic  conditions  at  the  wall. 
The  effect  on  the  coefficient  of  local  skin  friction 
of  varying  the  wall  temperature  was  not  included 
in  the  scope  of  the  present  series  of  tests.  Since 


this  effect  can  be  large  under  certain  conditions 
(ref.  19),  the  data  as  presented,  should  not  be  used 
directly  to  estimate  the  coefficient  of  local  skin 
friction  in  a  nonadiabatic  boundary  layer. 

EXPERIMENTAL  DETERMINATION  OF  THE  EQUIV.4LENCE  OF 

AIR  AND  HELIUM  IN  THE  TURBULENT  BOUNDARY  LAYER 

The  purpose  of  the  first  phase  of  the  testing 
program  was  to  evaluate  the  equivalence  of  air 
and  helium  in  the  turbulent  boundary  layer. 
This  was  clearly  necessary  prior  to  setting  up  an 
extensive  testing  program  using  helium  at  the 
higher  ^^ach  numbers.  The  evaluation  of  equiv¬ 
alence  was  made  from  comparisons  of  skin  fric¬ 
tion  and  profiles  for  boundary  layers  in  air  and 
helium.  The  air  boundary  layer  was  at  a  free- 
stream  Mach  number  of  approximately  4.2, 
wliereas  the  helium  boundary  layer  was  at  an 
actual  free-stream  Mach  number  of  approximately 
3.25,  or  at  an  equivalent  air  Mach  number, 
Ma,  of  approximately  4.2.  Since  tliis  comparison 
is  between  helium  and  air  at  wind-tunnel  tem¬ 
peratures  (Ti  —  fiOO®  R),  there  is  the  previously 
discussed  mismatch  in  m(Z')  between  the  two  gases 
in  the  boundarv  lavers  hit  the  wall,  «1.5 

^  '  -  ^"Alr 

Figure  16  shows  the  comparison  of  skin-friction 
results  obtained.  In  figure  16(a)  the  comparison 
is  made  between  skin  friction  for  air  and  helium 
boundarv  layers  as  a  function  of  the  Reynolds 
number  based  on  free-stream  viscosity  and  den¬ 
sity,  Rej.  (which  is  the  same  as  Rfi  at  the  skin- 
friction  element).  It  is  seen  that  the  local  skin- 
friction  coefficients  for  air  and  helium  fall  on 
essentially  one  curve.  In  figure  16(b),  where  the 
same  comparison  is  made  as  a  function  of  the 
Reynolds  number,  based  on  wall  viscosity 

and  density,  the  helium  points  fall  a  few  pei’cent 
above  the  air  points,  but  the  agreement  is  still 
reasonably  good.  Although  the  reasons  for  it 
are  not  known,  it  seems  clear  tluit  the  skin-fi*iction 
data  correlate  slightly  better  when  Reynolds 
number  is  based  on  free  stream  rather  than  on 
wall  properties.  (Consequently,  the  free-stream 
Reynolds  number  is  used  in  the  subsequent 
figures. 

Boundary -layer  profiles  for  both  air  and  helium 
at  an  equivalent  air  Mach  number  of  4.2  and  at 
equal  values  of  the  free-stream  Reynolds  number 
are  compared  in  figure  17.  The  profiles  were 
determined  from  surveys  taken  near  the  skin- 


Coetticient  or  local  sKin-rncfion,  Cf  ^  Coefftcienf  of  local  skin  friction, 
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Figure  16. — Equivalence  of  air  and  helium  turbulent  skin 
friction. 


friction  element.  In  some  cases  the  profiles  were 
determined  from  surveys  made  simultaneously  on 
either  side  of  the  element.  Data  from  these  sur¬ 
veys  are  shown  plotted  in  figure  17  (a)  and  (b)  and 
are  labeled  north  probe  and  south  probe.  An  aver¬ 
age  obtained  between  the  probes  is  given  in  figure 
17(c).  It  was  reasonable  to  use  this  average 
since  the  friction  velocity,  Ur,  was  calculated  from 
the  shear  sensed  by  the  skin-friction  element 
(which  was  not  at  the  exact  location  of  either 
probe) . 

The  data  in  figure  17  show  that  the  air  and 
helium  profiles  are  almost  identical.  In  parts  (a) 
and  (b)  of  the  figure,  the  small  differences  appear 
primarily  in  the  portions  of  the  profiles  near  the 
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wall.  This  is  not  surprising  since  with  equal 
free-stream  Reynolds  numbers  there  is  a  mismatch 
near  the  wall  in  the  dimensionless  viscosities  of 
air  and  helium.  In  part  (c)  of  the  figure  the 
differences  are  greatest  near  the  edge  of  the 
boundary  la3^er.  This  is  a  consistent  result  with 
that  of  parts  (a)  and  (b)  because  in  part  (c),  y* 
has  been  normalized  on  the  basis  of  a  wall  quantity, 
Ur,  and  differences  in  the  profiles  should  show  up 
most  prominently  away  from  the  wall.  Although 
the  air  and  helium  profiles  are  very  nearly  the 
same,  the  calculated  values  of  the  momentum 
thickness,  6,  do  not  correlate  perfectly  as  shown  by 
the  tabulation  in  figure  17.  This  variation  in  B  is 
of  the  same  order  as  the  traverse  variation  in  B  ob¬ 
tained  between  the  north  and  south  probes  and  is 
considered  to  be  equivalent  to  a  probable  error 
in  the  evaluation  of  Re^molds  number,  as  noted 
previously. 

Despite  the  calculated  differences  in  B  values, 
the  data  presented  in  figures  16  and  17  are  consid¬ 
ered  to  form  an  experimental  verification  of  the 
equivalence  of  air  and  helium  in  the  constant 
pressure  turbulent  boundary  layer.  The  good 
agreement  between  Of  values,  in  spite  of  the 
mismatch  of  m(2"),  indicates  that  the  turbulent 
skin-friction  coefficient  is  not  overly  sensitive  to 
variations  in  /u(T).  This  further  implies  that 
values  of  turbulent  Of  in  air  are  not  particularly 
sensitive  to  the  temperature  levels  of  the  air  and 
that  wind-tunnel  data  taken  with  air  should  be 
practically  valid  for  air  in  flight  (with  continuum 
conditions).  At  the  higher  Mach  numbers,  wind- 
tunnel  data  taken  with  helium  should  be  valid 
for  air  in  flight  since  the  ?(?)  match  between 
helium  and  flight  air  is  better  than  the  match 
between  helium  and  wind-tunnel  air. 

SKIN-FRICTION  RESULTS 

Using  the  concept  of  the  equivalence  of  air  and 
helium  in  the  turbulent  boundary  layer,  skin 
friction  data  were  taken  with  air  up  to  an  equiva¬ 
lent  air  Mach  number  of  4.2  and  with  helium  at 
the  higher  equivalent  air  Mach  numbers.  The 
results,  presented  in  several  ways,  are  summarized 
in  figures  18  throusrh  21,  inclusive.  Each  curve 
is  labeled  with  its  appropriate  equivalent  air 
Mach  number,  Ma. 

As  previously  discussed,  the  coefficient  of  local 
skin  friction  should  be  determined  as  a  function  of 
the  Reynolds  number  measured  from  the  virtual 
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origin,  Ecx  (with  the  other  parameters  held  fixed). 
However,  it  is  also  of  interest  to  observe  the  skin- 
friction  data  as  it  was  obtained,  with  the  Keynolds 
number  length  measured  from  a  fixed  point  in  the 
wind  tunnel.  These  data  are  presented  in  figure 
18  with  the  Reynolds  number,  measured  from 
the  nozzle  throat.  (For  consistency,  subsonic 
Reynolds  numbers  were  also  measured  from  this 
same  location  for  this  figure.)  Figure  18  shows 
the  type  of  data  obtained  when  the  Reynolds  num¬ 
ber  is  not  corrected  for  the  location  of  the  virtual 
origin,  and  it  also  shows  the  transition  characteris¬ 
tics  of  the  wind  tunnel.  The  solid  data  points  to 
the  left  of  the  peaks  show  transition  in  process  but 
not  completed;  these  points  do  not  represent  tur¬ 
bulent  data.  The  solid  points  to  the  right  of  the 
peaks  represent  turbulent  data,  but  with  transi¬ 
tion  less  than  8  inches  ahead  of  the  element  center. 
The  solid  points  are  not  considered  representative 
of  fully  turbulent  flows  and  are  not  shown  on  any 


other  figures.  As  mentioned  previously  the  flagged 
points  represent  turbulent  data  with  transition 
between  8  and  14  inches  ahead  of  the  element  cen¬ 
ter  (i?ez,— /?ex>0.40  Bej)  and  are  of  lesser  accu¬ 
racy  than  the  high  Reynolds  number  data  with 
transition  farther  upstream  (unflagged  points); 
these  points  are  also  flagged  on  all  subsequent 
figures.  The  curves  of  figure  18  are  not  considered 
to  represent  universal  relationships  between  local 
skin  friction  and  Reynolds  number  since  the  flow 
was  not  turbulent  over  the  entire  length  on  which 
Reynolds  number  was  based. 

With  the  Reynolds  number  measured  from  the 
virtual  origin  of  the  turbulent  boundary  layer,  the 
relationship  between  Of  and  Re^  is  the  universal 
one.  The  data  are  presented  in  this  form  in 
figure  19.  It  is  seen  that  the  flagged  points  tend  to 
become  somewhat  erratic  at  the  lowest  Reynolds 
numbers  (which  are  difficult  to  determine  accu¬ 
rately  as  previously  discussed) .  It  is  thought  that 
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Eicurb  19.— Varialioii  of  flat-plate  turbulent  skin  friction  with  Reynolds  number  based  on  lenp:{h  from  virtual  origin. 

wliich  faired  curves  representing  tlie  present  data 
are  compared  witli  faired  curves  representing  data 
taken  by  Scliiiltz-GriinoAv  (ref.  30),  Smitli  and 
Walker  (ref.  4),  and  Coles  (ref.  29).  It  is  noted 
tliat  the  comparison  is  reasonably  good  in  all  cases. 
At  M„=0.20,  tlic  present  data  and  the  subsonic 
data  of  Scliultz-Gnmow  and  Smitli  and  Walker 
fall  on  almost  the  same  curve. 

In  perhaps  the  most  convenient  form,  the  pres¬ 
ent  data  are  summarized  in  figure  21  showing  tlie 
ratio  of  the  local  skin-friction  coefficient  to  the  in¬ 
compressible  local  skin-friction  coefficient  (Kar- 
man-Schoenherr)  at  tlie  same  Reynolds  number, 
CfjCf.,  as  a  function  of  Reynolds  number, 
Comparison  of  the  levels  of  the  curves  shows  the 
large  effect  of  Alach  number  on  CfjCf..  Any 
Reynolds  number  effect  on  CfjCf.  for  a  given  Mach 
number  can  be  seen  to  be  small.  In  examining  the 
unflagged  points,  it  is  difficult  to  discern  any  strong 
Reynolds  number  trend  within  tlie  scatter  of  the 
experimental  data.  It  appears  that  for  all  prac¬ 
tical  purposes,  CfjCf.  can  be  regarded  as  inde¬ 
pendent  of  Reynolds  number.  This  is  the  same 
result  as  tliat  previously  obtained  on  bodies  of 
revolution  by  Chapman  and  Kester  (ref.  18)  over 
a  more  limited  range  of  test  conditions. 


the  unflaggcd  points  provide  accurate  values  of 
the  coefficient  of  local  skin  friction  for  the  Mach 
numbers  and  Reynolds  numbers  at  which  the  data 
were  obtained.  The  subsonic  data  arc  in  close 
agreement  with  the  Karin an-Schoenherr  curve  for 
incompressible  flow  which  has  been  included  here 
for  comparison  purposes. 

The  present  data  have  been  compared  with  that 
taken  by  other  investigators.  The  ranges  of  Mach 
numbers  and  Reynolds  numbers  over  which  com¬ 
parison  was  possible  arc  shown  in  figure  20  in 


Reynolds  number,  RBji,  millions 


Fictjre  20. — Comparison  of  direct  force  moasiiremonts  of 
turbulent  skin  friction. 
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COMPARISON  OF  EXPERIMENTAL  RESULTS  WITH  METHODS 
OF  CALCULATION 

The  experimental  values  obtained  for  the  eoeffi- 
eient  of  local  skin  friction  were  compared  witli  tlie 
results  of  four  metJjods  of  calculating  local  skin 
friction  in  the  compressible  turbulent  boundary 
layer:  tlie  T'  metliod  (ref.  19),  the  metJiod  of  Walz 
(ref.  31),  the  Van  Driest  method  using  the  Prandtl 
mixing  length  (ref.  21),  and  the  Van  Driest  method 
using  the  Von  Karman  mixing  length  (ref.  32). 
Tlie  four  methods  are  not  necessarily  considered 
‘‘best^'  by  tlie  authors,  but  were  selected  as  being 
representative  examples  from  the  large  numbrr  of 
methods  proposed  for  calculating  skin  friction  in 
the  compressible  turbulent  boundaiy  layer.  The 

-  T' 

- Walz 

- Van  Driest  (Prandtl  mixing  length) 

- Von  Driest  (Karman  mixing  length) 


Flagged  symbols,  Re^/Rei_<0.^0 


—  - 

— 

r" 

ir: 

— 

- 

—  - 

rr-  O' 

o-  ■ 

— 

(a) 

>nP<>-CK 

.5 


A 


— 

_ _ 

• — 

— ■ 

— 

(d) 

_ 

■3^ 

AS-- 

bz: 

— 

— 

rz 

I  2  3  4  6  8  10  20  30  40  60  80  IOC 

Reynolds  number,  Rex,  millions 

(a)  i\/„  =  2.95 

(b)  A/„=4.2 

(c)  A/, -6.7 

(d)  A/„  =  9.9 

Figure  22. — Calculated  C/IC/i  coniparcd  with  experi¬ 
mental  CfICfi. 


comparisons  are  shown  in  figure  22.  In  calculatin 
tlie  ordinate,  O/IC/.,  the  value  in  all  cases  wa 
taken  as  tlie  Karman-Schoeidierr  value. 

Jt  is  noted  that  of  the  four  methods  the  T'  (in¬ 
termediate  enthalpy)  method  appears  to  give  the 
best  fit  to  the  experimental  data.  The  T'  method 
uses  a  means  reference  temperature  to  evaluate 
the  physical  piopertii's  of  the  gas.  It  was  oi'igi- 
nally  used  by  Rubesin  and  eJohnson  (ref.  33)  for 
laminar  flows  and  was  adapted  to  turbulent  flows 
by  Eckert,  Sommei’  and  Short,  and  others.  The 
constants  used  in  the  cahnilations  are  those  deter¬ 
mined  by  Sommer  and  Short  (ref.  19).  (Constants 
obtained  by  other  investigators  do  not  change  the 
results  much.)  It  is  seen  that  th(‘  experimental 
data  and  the  T'  calculations  agree  reasonably  well 
although  the  trends  with  Reynolds  number  are 
not  in  complete  agreement.  Tlie  results  of  the  T' 
calculations  are  somewhat  aflected  by  the  viscosity 
temperature  relationship  used  for  the  gas  in  the 
boundary  layer.  This  means  that  the  results  for 
air  and  helium  will  differ  slightly.  The  results  for 
air  are  also  affected  by  the  temperature  levels  of 
th('  air  {Ti  or  To^)  siiu'e  the  dinu'usionh'ss  viscosity 
temperature  rc'lationship  for  air  depends  on  the 
temperature  levcd  used.  The  T'  curves  in  figure  22 
were  calculated  for  air  at  wind-tunnel  tempera¬ 
tures  to  compare  with  data  taken  at  equivalent  air 
Mach  numbers  of  2.95  and  4.2;  the  T'  curves 
were  cadculated  for  hc'lium  to  compare  with  data 
taken  at  equivalent  air  Mach  numlicrs  of  6.7  and 
9.9.  TJu'se  calculations  corresjiond  to  the  actual 
running  conditions  in  the  bound  ary -laycT  channed. 
TJie  effects  on  the  T'  results  of  different  tempera¬ 
ture  levcds  (with  air)  and  of  changing  from  air  to 
hedium  are  shown  in  figure  23.  On  a  percentage 
(but  not  absolute)  basis,  the  spread  of  these  Cy/Cy. 
curve's  increases  with  incTease  in  Macdi  number. 

C^alculations  of  the  coefficient  of  the  local  skin 
friction  by  the  method  of  Walz  are  also  shown  in 
figure  22.  This  is  a,  semic'mpirical  intc'gral  method 
that  makes  use  of  results  obtained  from  calcaila- 
tions  of  incompressible  flows.  It  is  sc'en  in  the 
figure  that  some  values  of  calculated  by 

this  method  are  higher  and  some  are  lower  than 
the  experimental  points,  but,  in  all  cases,  th(\v 
increase  with  increasing  Reynolds  number  at  a 
given  Mach  number.  This  trend  is  not  strong, 
but  the  trend  of  the  curves  does  not  fit  the  e.x- 
perimental  data  closely.  By  a  visual  extrapola- 
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- Air  in  wind  tunnel,  Tf  =  600  ®R 

-  Air  in  flight,  T/ -(l+0.2/l^^)(400),  ®R 

- Helium 


Reynolds  number,  millions 

Fictre  23. — Effect  of  temperature  levels  and  gas  proper¬ 
ties  on  skin  friction  calculated  by  the  T'  method. 

tion,  it  can  be  seen  that  these  calculated  curves 
will  cross  the  curves  of  experimental  data  at  some 
Reynolds  number. 

As  seen  in  figure  22,  the  values  of  Cf/Cf.  calcu¬ 
lated  by  the  Van  Driest  method  with  Prandtl 
mixing  length  run  consistently  higher  than  the 
experimental  data.  Also  the  CfjCf.  values  fall 
off  fairly  rapidly  with  increasing  Reynolds  num¬ 
ber,  whereas  the  experimental  data  do  not  indicate 
this  effect.  Values  calculated  by  the  Van  Driest 
method  with  Von  Karmen  mixing  length  are 
considerably  closer  to  the  experimental  data 
(although  still  a  little  high),  and  the  slopes  of  the 
curves  correspond  more  closely  to  the  trend  of  the 
data. 

SUBSONIC  DATA  AND  POSSIBLE  ROUGHNESS  EFFECTS 

Subsonic  data  were  taken  principally  to  deter¬ 
mine  the  accuracy  and  validity  of  the  experimental 
techniques  used,  since  these  data  could  be  com¬ 
pared  with  the  well-known  properties  of  the  low- 
speed  turbulent  boundary  layer.  Specifically,  it 
was  desired  to  find  out  whether  the  supersonic 
data  could  possibly  contain  any  roughness  effects. 
It  has  been  shown  that  high  Mach  number 
boundary  layers  are  less  sensitive  to  roughness 
than  are  subsonic  boundary  layers  (ref.  34).  The 
subsonic  and  supersonic  data  were  taken  with 
the  same  top  plate  and  skin-friction  element,  so 
that  the  roughness,  if  any,  was  invariant.  Any 
roughness  effects,  then,  would  appear  most 
strongly  in  the  subsonic  data.  As  is  well  known, 
when  appreciable  roughness  exists,  the  subsonic 


skin-friction  curve  should,  with  increasing  Reyn¬ 
olds  number  for  a  fixed  length,  diverge  upward 
from  the  curve  for  a  smooth  surface  and  eventually 
become  horizontal.  As  seen  in  figure  19,  the 
subsonic  data  exhibit  no  such  roughness  effects. 
It  is  therefore  believed  that  the  supersonic  data 
are  free  of  roughness  effects. 

BOUNDARY-LAYER  PROFILES 

Although  the  main  purpose  of  this  research  was 
to  obtain  and  report  skin-friction  measurements, 
it  was  clearly  of  interest  also  to  secure  some 
boundary-layer  profiles  which  would  at  least 
provide  a  qualitative  “look’’  at  the  boundary 
layers  being  studied.  Typical  profiles  are  shown 
in  figures  24,  25,  and  26,  the  coordinates  in  figure 
24  being  M/M^  against  y/O,  the  coordinates  in 
figure  25  being  u/U^  against  y/d  for  the  same 
profiles,  while  in  figure  26  the  same  profiles  again 
were  plotted  with  u*  against  y*. 

The  MjM^  against  yjd  plots  in  figure  24  repre¬ 
sent  data  directly  reduced  from  pitot  traverses. 
It  is  seen  that  all  of  the  M/M„  profiles  are  some¬ 
what  similar  in  shape;  however,  some  variation  in 
shape  with  change  in  Mach  number  can  be  seen. 
A  decrease  in  Bjh  with  increase  in  Mach  number  is 
readily  recognized.  Also  it  is  seen  that  MjM^  at 
the  edge  of  the  sublayer  decreases  with  an  increase 
in  Mach  number.  In  physical  dimensions  the 
sublayer  becomes  thicker  with  an  increase  in 
Mach  number  although  this  is  not  readily  seen 
from  the  normalized  coordinates  used  in  the  plots. 
No  corrections  were  made  to  the  pitot  readings  to 
account  for  the  fluctuating  (turbulent)  component 
of  velocity  since  it  was  felt  that  these  corrections 
would  be  small  and  would  be  less  than  the  ex¬ 
perimental  scatter  of  the  data.  If  these  correc¬ 
tions  had  been  made,  the  Mach  numbers  near 
the  wall  would  have  been  reduced  by  a  few  percent 
since  the  pitot  probe  senses  values  higher  than  the 
time  mean  in  turbulent  flow.  Also,  as  discussed 
below,  it  is  believed  that  the  largest  error  by  far, 
is  that  of  the  distortion  of  the  boundary  layer  due 
due  to  the  presence  of  the  probe.  As  observed 
previously  in  discussion  of  figure  17,  a  transverse 
variation  in  the  value  of  B  of  approximately  11 
percent  can  also  be  seen  in  figure  24  (a)  which  shows 
surveys  taken  on  either  side  of  the  skin-friction 
element.  The  profiles  again  fall  together  when 
plotted  on  a  yjB  basis. 

The  profiles  in  figure  25  convey  essentially  the 
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same  information,  in  somcwliat  different  form,  as 
the  profiles  in  figure  24.  The  velocity,  u,  was 
calculated  from  the  pitot  data  by  assuming  a 
constant  total  temperature  (settling  chamber 
temperature)  through  the  boundary  layer.  This 
assumption,  slightly  in  error,  gives  calculated 
velocities  as  much  as  approximately  5  percent 
too  high  near  the  wall  (assuming  the  pitot  Mach 
number  is  correct).  With  profiles  plotted  on 
a  dimensionless  velocity  basis,  the  points  are 
more  spread  out  near  the  wall  (in  the  sublayer) 
than  when  plotted  on  an  MjM^  basis,  as  most  of 
the  velocity  gradient  occurs  in  the  wall  portion 
of  the  boimdaiy  layer. 

Also  sliown  on  each  profile  in  figures  24  and  25 
is  a  straight  line  from  the  origin  with  slope  deter¬ 
mined  by  the  directly  measured  coefficient  of  local 
skin  friction.  (Where  profiles  from  north  and 
south  profiles  are  superposed,  average  values  of  0 
were  used  in  calculating  the  slope  of  the  straight 
line.)  It  is  noted  that  with  increasing  Mach  num¬ 
bers  there  is  an  increasing  discrepancy  between 
the  slope  deduced  from  skin-friction  measurements 
and  the  apparent  slope  of  tlie  measured  profiles. 
This  difference  is  presumably  not  due  to  any 
transverse  nonuniformities  in  the  flow  as  this 
difference  was  sensed  by  probes  on  either  side  of 
the  skin-friction  element.  It  is  believed  that  tlie 
direct  measurements  of  skin  friction  provide  tlie 
more  accurate  determination  and  that  the  ob¬ 
served  difference  is  largely  duo  to  tlie  distortion  of 
the  sublayer  by  the  presence  of  the  probe  in  the 
sublayer.  This  efl'ect  may  be  three-dimensional, 
with  the  slower  streamlines  near  the  wall  being 
displaced  sidewards  and  the  faster  streamlines 
(located  away  from  tlie  wall)  being  displaced  to¬ 
ward  the  wall.  This  could  account  for  a  probe 
near  the  wall  sensing  higher  velocities  than  would 
normally  be  expected.  Similar  distortions  in 
turbulent  incompressible  flows  have  been  observed 
by  a  number  of  investigators.  G.  I.  Taylor  re¬ 
ported  on  experiments  with  a  Stanton  tube  show¬ 
ing  a  Reynolds  number  effect  on  the  distortion 
(ref.  35).  Von  Doenhoff  showed  a  similar  Rey¬ 
nolds  number  effect  with  a  total  pressure  tube  in 
contact  with  a  surface  (ref.  36).  There  are  very 
few  data  for  probes  located  away  from  the  wall. 
Young  and  Maas  (ref.  37)  made  pitot  measure¬ 
ments  in  a  wake,  but  their  conclusion  that  the 
streamline  displacement  is  approximately  one- 
fourth  of  the  probe  diameter  seems  too  small  for 


the  turbulent  boundai'v  layer.  As  pointed  out  by 
Von  Doenhoff,  the  Young  and  Maas  displacement 
is  not  large  enough  to  be  consistent  with  values 
measured  at  the  wall  (and  it  should  also  have  a 
Reynolds  number  dependence  as  determined  from 
measurements  at  the  wall). 

The  apparent  distortion  of  the  boundary  laycj’ 
near  the  wall  observed  by  the  authors  is  of  the 
same  order  of  magnitude  as  tliat  reported  by 
Taylor  and  Von  Doenhoff.  In  noting  the  dis¬ 
agreement  between  the  measured  profile  slope 
near  the  wall  and  the  slope  deduced  from  the  skin 
friction,  it  is  felt  that  caution  should  be  exer¬ 
cised  in  placing  a  high  degree  ot  reliance  on  pitot 
data  in  compressible  turbulent  flows  near  a  wall. 

In  figure  26  the  same  data  arc  plotted  showing 
u'^  —  ulvr  against  y'^  =  yuTlvjr.  The  north  and 
south  probes  are  not  distinguished  in  this  figure, 
and  in  the  cases  where  two  profiles  were  obtained 
simultaneously,  average  values  of  n  and  y  were 
used  in  calculating  'W*  and  ?/*.  Also,  because  of 
overcrowding,  points  near  the  outer  edge  of  the 
boundary  layer  were  thinned  out  wliere  necessary 
for  this  figure.  For  purposes  of  comparison,  tlie 
curves,  =  and  v,*  =  5.75  logioy*  +  5.10  are 
shown  in  figure  26.  The  constants  in  tlie  latter 
equation  (for  incompressible  flows)  are  obtaitied 
from  reference  38.  In  figure  26(g)  are  plotted  one 
profile  foi*  eacli  ^^ach  number,  all  at  a]}proximately 
the  same  Reynolds  number.  It  is  seen  that  all  of 
these  profiles  have  approximately  the  same  shape. 
With  the  small  amount  of  data  obtained  in  the 
sublayer  and  at  the  edge  of  tlie  sublayer,  it  does 
not  appear  to  be  possible  to  distinguish  a  definite 
trend  with  Mach  number  of  the  variation  oi  v* 
and  ?/*  at  the  edge  of  the  sublayer. 

If  the  profiles  in  figures  24  and  25  were  to  be 
arbitrarily  shifted  to  correspond  to  the  measured 
skin-friction  values,  a  similar  correction  would  be 
made  in  the  1/*  vs.  ?/*  profiles  in  figure  26.  This 
would  correspond  to  the  so-called  Nikuradse  shift 
and  would  be  of  approximately  the  same  order  of 
magnitude  (ref.  39).  As  reported  in  reference  39 
Nikuradse  arbitrarily  shifted  liis  -?/*  vs.  ?/*  profiles 
outward  by  a  constant  amount  so  that  the  points 
nearest  the  wall  would  be  consistent  with  the  skin 
friction.  In  the  absence  of  data  giving  accurate 
values  of  probe-induced  distortion  throughout  the 
boundary  layer,  it  would  appear  that  a  shift  of 
the  entire  profile  is  not  waiTant(‘d.  However,  it 
should  be  kept  in  mind  that  the  portion  of  tlie 
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(e)  Ma=6.7;  helium;  Re,=7.2XW 


(f)  M„=9.9;  helium;  fie^=7.3X106  (C/ flagged) 


Figure  26. — Continued 
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(g)  Mach  numbers  indicated. 
Figuke  26. — Concluded 


Turbulent  profile -helium.  (/?e.,.//T’e^<0.6o) 

/?fi-^  =  l8.7xl06 

O  J  Transitional  profile-helium 

l/P/^  =  10.6  18.0X10® 

- Direct  skin-friction  measurement 

- (Laminar  profile-calculated  (Hantzsche  and  Wendt) 

}Ma  =  \0,6  18.0X10® 


0  .1  .2  .3  .4  .5  .6  .7  .8 


/,  inches 

Figure  27.- — Comparison  of  laminar,  turbulent,  and  tran¬ 
sitional  profiles. 

profile  near  the  wall  becomes  somcwliat  distorted 
by  the  presence  of  a  probe,  and  the  data  points 
near  the  wall  are  probably  not  reliable  while  tliose 


away  from  the  wall  probably  are.  Shifting  tlie 
entire  profile  by  a  constant  amount  would  s(‘arcdy 
be  noticed  for  points  far  removed  from  the  wall, 
blit  it  would,  of  course,  aflV'ct  the  points  near  the 
wall. 

As  a  matter  of  interest,  figure  27  has  been  in¬ 
cluded  to  show  a  comparison  between  a  turliulent 
profile,  a  transitional  profile,  and  a  calculated 
laminar  profile  (ref,  40),  all  at  approximately  the 
same  Reynolds  number  from  the  nozzle  tlnoat, 
/iVi,  and  all  at  equivalent  air  Mach  numbers  near 
10.  It  is  of  interest  to  note  that  the  slope  of  the 
straight  line  deduced  from  skin-friction  measure¬ 
ments  checks  reasonably  well  with  the  wall  slope 
of  the  transitional  profile.  As  lias  been  pointed 
out,  these  slopes  did  not  check  each  other  with 
the  high  Macli  number  tiu’bulent  profih's  that  were 
obtained,  including  the  one  shown  in  figure  27,  It 
would  appear  that  the  probe  produces  less  distor¬ 
tion  near  the  wall  for  transitional  boundary-layer 
profiles  than  for  turbulent  profiles.  The  explana¬ 
tion  for  this  is  not  known. 

CONCLUSIONS 

Skin  friction  was  measured  in  constant-pressure, 
adiabatic,  turbulent  boundary  layers  up  to  a 
maximum  Reynolds  number  of  100X10*^  (at 
Ma=4.2)  and  up  to  a  maximum  equivalent  air 
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Mach  number  of  9.9.  Air  and  helium  were  used 
as  test  media.  From  the  data  obtained,  the  fol¬ 
lowing  conclusions  have  been  drawn : 

1.  It  is  possible  to  set  up  equivalent,  constant- 

pressure,.  turbulent  boundary  layers  using  different 
gases  as  test  media.  Measurements  of  skin  fric¬ 
tion  taken  with  helium  agree  closely  with  those 
taken  with  air  when  is  the  same  for 

the  two  gases  and  wall  conditions  are  adiabatic. 

2.  The  variation  of  the  coefficient  of  local  skin 
friction  as  calculated  by  the  T'  method  agrees 


particularly  well,  while  the  Van  Driest  method 
with  Von  K4rm^n  mixing  length  agrees  reasonably 
well  with  experimental  data  over  the  range  of  test 
conditions  investigated. 

3.  The  experimental  data  indicate  that  the 
ratio,  is  essentially  independent  of  Reynolds 

number. 

Ames  Research  Center 

National  Aeronautics  and  Space  Administration 
Moffett  Field,  Calif.,  June  20,  1960 


APPENDIX  A 

EXPERIMENTAL  DETERMINATION  OF  TEMPERATURE  RECOVERY  FACTOR 


Adiabatic  wall  conditions  in  the  turbulent 
boundary  layers  were  established  experimentally 
without  actually  knowing  the  temperature  recov¬ 
ery  factor.  However,  it  is  of  interest  to  calculate 
the  recovery  factors  as  follows: 


Vr  = 


(Al) 


In  more  convenient  form: 


Vr  = 


2 


(A2) 


Experimentally,  tlie  adiabatic  wall  temperature, 
Tawy  was  taken  as  the  existing  wall  temperature 
prior  to  the  run.  The  Mach  number  was  known 
prior  to  the  run  from  previous  measurements  with 
the  given  nozzle  and  test  section  geometry.  The 
total  temperature,  Tt,  was  adjusted  before  the 
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run  by  setting  the  desii’ed  temperature  in  the 
heat  exchanger.  In  operation  the  heat  exchanger 
held  the  Tt  constant  within  several  degrees  until 
tlie  heat  capacity  of  the  exchanger  was  used  up, 
at  which  time  tlie  7'^  fell  off  rapidly,  and  the  run 
was  then  terminated.  The  adiabatic  condition 
was  determined  by  the  constancy  of  the  wall 
temperature  during  the  run.  A  2°  F  change  in 
wall  temperature  during  a  run  was  considered 
allowable  for  adiabatic  operation. 

Values  of  Tt  were  initially  preset  by  calculation 
from  expected  values  of  Adjust  in  ents  of  Tt 

for  subsequent  runs  were  based  on  experience, 
with  constancy  of  wall  temperature  as  the  cri¬ 
terion.  Recovery  factors  were  then  calculated 
from  the  acceptable  values  of  Tt. 

Recovery  factors  were  0.87  for  air  at  Ma  values 
of  2.95  and  4.2.  For  helium  the  recoviuy  factors 
were  0.80,  0.87,  and  0.88  for  values  of  4.2, 
6.7,  and  9.9,  respectively. 


APPENDIX  B 

ELECTRICAL  CIRCUITRY  OF  THE  SKIN-FRICTION  ELEMENT  AND  SERVO 


The  electrical  circuitry  of  the  skin-friction 
element  and  servo  was  designed  and  the  equipment 
constructed  under  the  supervision  of  William  J. 
Kerwin  of  the  Ames  Research  Center.  The 
skin-friction  element  transducer  consists  of  two 
variable  inductors  whose  inductances  are  changed 
in  opposite  directions  by  the  motion  of  a  mag¬ 
netic  armature  (magnetic  return  piece)  on  the 
disk  (see  figs.  8  and  9).  The  armature  varies  the 
size  of  a  small  air  gap  and  thus  changes  the  reluc¬ 
tance  of  the  magnetic  path. 

Figure  10  shows  the  double  balance  a-c  bridge 
and  associated  circuitry  used.  The  relations 
necessary  for  balance  are: 

Tji  Ro“1“0!R 

l;“Ro+(i-«)r  ^  ^ 


where 

go  gap  spacing  with  the  element  disk  in  the 
neutral  position 

g  variation  in  gap  spacing  from  the  neutral 
position 

Ki  constant  of  proportionality 


Equation  (Bl)  then  reduces  to: 

go+g  1^0+  (1 — «)  ^ 


Solving  for  a  one  obtains 


1 

^”2  V  2goR  ; 


g 


— KgAg 


(B5) 


(B6) 

(B7) 


ri  Ro”|“aR 

r2“Ro+(l-a)R  ^  ^ 

where 

Li  variable  inductance  of  skin-friction  element 
coil  No.  1 

L2  variable  inductance  of  skin-friction  element 
coil  No.  2 

Ro  fixed  resistances  in  the  bridge 
R  total  resistance  of  potentiometer  No.  1  (con¬ 
stant) 

a  fractional  part  of  R 

ri  total  resistance  (variable)  in  skin-friction 

element  coil  No.  1  (fixed)  and  partial  resist¬ 
ance  of  potentiometer  No.  2  (variable) 
r2  total  resistance  (variable)  in  skin-friction 

element  coil  No.  2  (fixed)  and  partial  re¬ 
sistance  of  potentiometer  No.  2  (variable) 
Since  the  inductances  are  approximately  inversely 
proportional  to  the  air  gap  spacing,  one  can  write : 


(B3) 


L2- 


^_i_ 

go-g 


(B4) 


where 

R  +  2R0 
2goR 

Thus  the  reading  of  a  counter  connected  to  poten¬ 
tiometer  No.  1  is  directly  proportional  to  the 
change  in  the  gap  spacing  in  the  skin-friction 
element. 

The  bridge  balance  is  accomplished  by  an  auto¬ 
matic  dual  balance  servo.  One  servomotor  is 
connected  to  potentiometer  No.  1  and  is  driven  so 
as  to  null  the  component  of  bridge  output  voltage 
in  quadrature  with  the  applied  bridge  voltage, 
this  component  being  introduced  primarily  by 
inductive  unbalance  of  the  bridge  (eq.  (Bl)).  A 
second  servomotor  is  connected  to  potentiometer 
No.  2  and  is  driven  so  as  to  null  the  component 
of  bridge  output  voltage  which  is  in  phase  with 
the  applied  bridge  voltage,  this  component  being 
introduced  primarily  by  resistance  unbalance  of 
the  bridge  (eq.  (B2)).  A  3000  count  dial  con¬ 
nected  to  the  shaft  of  potentiometer  No.  1  indicates 
the  displacement  of  the  disk  in  the  skin-friction 
element. 
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